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^SJ ' Abstract 

^ ' Design problems in industrial engineering often involve a large number of design variables 

with multiple objectives, under complex nonlinear constraints. The algorithms for multiobjec- 
tive problems can be significantly different from the methods for single objective optimization. 
To find the Pareto front and non-dominated set for a nonlinear multiobjective optimization 
problem may require significant computing effort, even for seemingly simple problems. Meta- 
heuristic algorithms start to show their advantages in dealing with multiobjective optimization. 
In this paper, we extend the recently developed firefly algorithm to solve multiobjective opti- 
, mization problems. We validate the proposed approach using a selected subset of test functions 

' and then apply it to solve design optimization benchmarks. We will discuss our results and 

• ' provide topics for further research. 
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1 Introduction 



Design optimization in engineering and industry often concerns multiple design objectives under 
complex, highly nonlinear constraints. Different objectives often conflict each other, and sometimes, 
truly optimal solutions do not exist, and some compromises and approximations are often needed 
[1, 2, 3]. Further to this complexity, a design problem is subjected to various design constraints, 
limited by design codes or standards, material properties and the optimal utility of available resources 
and costs [2, 4]. Even for global optimization problems with a single objective, if the design functions 
are highly nonlinear, global optimality is not easy to reach. Metaheuristic algorithms are very 
powerful in dealing with this kind of optimization, and there are many review articles and textbooks 
[5, 6, 7, 8, 9, 10, 11]. 

In contrast with single objective optimization, multiobjective problems are much more diflacult 
and complex [5, 12]. Firstly, no single unique solution is the best; instead, a set of non-dominated 
solutions should be found in order to get a good approximation to the true Pareto front. Secondly, 
even if an algorithm can find solution points on the Pareto front, there is no guarantee that multiple 
Pareto points will distribute along the front uniformly, often they do not. Thirdly, algorithms which 
work well for single objective optimization usually cannot directly work for multiobjective problems, 
unless under special circumstances such as combining multiobjectives into a single objective using 
some weighted sum method. Substantial modifications are needed to make algorithms for single 
objective optimization work. In addition to these difficulties, a further challenge is how to generate 
solutions with enough diversity so that new solutions can sample the search space eflSciently. 

Furthermore, real-world optimization problems always involve some degree of uncertainty or 
noise. For example, materials properties for a design product may vary significantly. An optimal 
design should be robust enough to allow such inhomogeneity, which provides a set of multiple 
feasible solution sets. Consequently, optimal solutions among the robust Pareto set can provide good 
options so that decision-makers or designers can choose to suit their needs. Despite these challenges, 
multiobjective optimization has many powerful algorithms with many successful applications [6, 13, 
14, 15, 43]. In addition, metaheuristic algorithms start to emerge as a major player for multiobjective 
global optimization, they often mimic the successful characteristics in Nature, especially biological 
systems [9, 10], while some algorithms are inspired by the beauty of music [16]. Many new algorithms 
are emerging with many important applications [8, 10, 13, 17, 18, 19, 20]. 

For example, multiobjective genetic algorithms are widely known [15, 21], while multiobjective 
differential evolution algorithms are also very powerful [22, 23]. In addition, multiobjective particle 
swarm optimizers are becoming increasingly popular [19]. As there are many algorithms, one of our 
motivations in the present study is to compare the performance of these algorithms for real-world 
application. 

Most metaheuristic algorithms are based on the so-called swarm intelligence. PSO is a good 
example, it mimics some characteristics of birds and fish swarms. Recently, a new metaheuristic 
search algorithm, called Firefly Algorithm (FA), has been developed by Yang [9, 11]. FA mimics 
some characteristics of tropic flrefly swarms and their flashing behaviour [10, 11]. A firefly tends to 
be attracted towards other fireflies with higher flash intensity. This algorithm is thus different from 
PSO and can have two advantages: local attractions and automatic regrouping. As light intensity 
decreases with distance, the attraction among flreflies can be local or global, depending on the 
absorbing coefficient, and thus all local modes as well as global modes will be visited. In addition, 
fireflies can also subdivide and thus regroup into a few subgroups due to neighboring attraction is 
stronger than long-distance attraction, thus it can be expected each subgroup will swarm around a 
local mode. This latter advantage makes it particularly suitable for multimodal global optimization 
problems [10, 11]. 

Preliminary studies show that it is very promising and could outperform existing algorithms 

such as particle swarm optimization (PSO). For example, a Firefly-LGB algorithm, based on firefly 
algorithm and Linde-Buzo-Gray (LGB) algorithms for vector quantization of digital image compres- 
sion, was developed by Horng and Jiang [24] , and their results suggested that Firefly-LGB is faster 
than other algorithms such as particle swarm optimization LBG (PSO-LBG) and honey-bee mating 
optimization LBG (HBMO-LBG). Apostolopoulos and Vlachos provided a detailed background and 
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analysis over a wide range of test problems [25], and they also solved multiobjective load dispatch 
problem using a weighted sum method by combining multiobjectives into a single objective, and 
their results are very promising. The preliminary successful results of firefly-based algorithms pro- 
vide another motivation for this paper. That is to see how this algorithm can be extended to solve 
multiobjective optimization problems. 

In this paper, we will extend FA to solve multiobjective problems and formulate a multiobjective 
firefly algorithm (MOFA). We will first validate it against a subset of multiobjective test functions. 
Then, we will apply it to solve design optimisation problems in engineering, including bi-objective 
beam design and a design of a disc brake. Finally, we will discuss the unique features of the proposed 
algorithm as well as topics for further studies. 

2 Multiobjective Firefly Algorithm 

In order to extend the firefly algorithm for single objective optimization to solve multiobjective 
problems, let us briefly review its basic version. 

2.1 The Basic Firefly Algorithm 

Firefly Algorithm was developed by Yang for continuous optimization [9, 10, 27], which was subse- 
quently applied into structural optimization [26] and image processing [24]. FA was based on the 
flashing patterns and behaviour of fireflies. In essence, FA uses the following three idealized rules: 
(1) Fireflies are unisex so that one firefly will be attracted to other fireflies regardless of their sex; (2) 
The attractiveness of a firefly is proportional to its brightness and they both decrease with distance. 
Thus for any two fiashing fireflies, the less brighter one will move towards the brighter one. If there 
is no brighter one than a particular firefly, it will move randomly; (3) The brightness of a firefly is 
determined by the landscape of the objective function. 

For a maximization problem, the brightness can simply be proportional to the value of the 
objective fimction. As both light intensity and attractiveness affect the movement of flreflies in the 
firefly algorithm, we have to define their variations. For simplicity, we can always assume that the 
attractiveness of a firefly is determined by its brightness which in turn is associated with the encoded 
objective function. In the simplest case for maximum optimization problems, the brightness / of a 
firefly at a particular location x can be chosen as I{x) oc ,f{x). 

However, the attractivcnicss 13 is relative, it shoiild be seen in the eyes of the beholder or judged 
by the other firefiies. Thus, it will vary with the distance between firefiy i and firefiy j. Therefore, 
we can now define the attractiveness ^ of a firefly by 

P = l3oe-'''\ (1) 

where /3o is the attractiveness at r = 0. In fact, equation (1) defines a characteristic distance 
F = over which the attractiveness changes significantly from /3o to (3oe~^- The distance 

between any two fireflies i and j at Xi and xj , respectively, is the Cartesian distance r jj = 1 1 x j — a; ^ 1 1 . 
It is worth pointing out that the distance r defined above is not limited to the Euclidean distance. 
In fact, any measure that can effectively characterize the quantities of interest in the optimization 
problem can be used as the 'distance' r. We can define other distance r in the n-dimensional 
hyperspace, depending on the type of problem of our interest. 

For any given two firefiies Xi and Xj, the movement of firefly i is attracted to another more 
attractive (brighter) firefly j is determined by 

xl+' ^x\+ /3oe-^'-- (4 - x'^) + at e\, (2) 

where the second term is due to the attraction. The third term is randomization with at being the 
randomization parameter, and e* is a vector of random numbers drawn from a Gaussian distribution 
or uniform distribution. The location of fireflies can be updated sequentially, by comparing and 
updating each pair of them in every iteration cycle. 
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For most implementations, we can take /3o = 1 and at = 0(1), though we found that it is better 
to use a time-dependent at so that randomness can be reduced gradually as iterations proceed. It 
is worth pointing out that (2) is a random walk biased towards the brighter fireflies. If /3o = 0, it 
becomes a simple random walk. Furthermore, the randomization term can easily be extended to 
other distributions such as Levy flights [28]. 

The parameter 7 now characterizes the variation of the attractiveness, and its value is crucially 
important in determining the speed of the convergence and how the FA algorithm behaves. In 
theory, 7 G [0,oo), but in practice, 7 = 0(1) is determined by the characteristic distance F of the 
system to be optimized. Thus, for most applications, it typically varies from 10~^ to 10^. 

To consider the scale variations of each problem, we now use rescaled, vectorized parameters 

a = 0.01L, 7 = 0.5/L2, (-g-, 

with L = {Ub — Lb) where Ub and Lb are the upper and lower bounds of a;, respectively. Here the 
factor 0.01 is to make sure the random walks is not too aggressive, and this value has been obtained 
by a parametric study. 

2.2 Multiobjective Firefly Algorithm 

For multiobjective optimization, one way is to combine all objectives into a single objective so that 
algorithms for single objective optimization can be used without much modifications. For example, 
FA can be used directly to solve multiobjective problems in this manner, and a detailed study was 
carried out by Apostolopoulos and Vlachos [25]. 

Another way is to extend the firefly algorithm to produce Pareto optimal front directly. By 
extending the basic ideas of FA, we can develop the following Multi-objective Firefly Algorithm 
(MOFA), which can be summarized as the pseudo code listed in Fig. 1. 



Define objective functions fi{x), fK{x) where x = {x\, Xd)^ 
Initialize a population of n fireflies a;, (i = 1, 2, n) 
while {t <MaxGeneration) 
for i,j = \:n (all n fireflies) 

Evaluate their approximations PFi and PFj to the Pareto front 

\H ^ j and when all the constraints are satisfied 
if PFj dominates PFi, 

Move firefiy i towards j using (2) 

Generate new ones if the moves do not satisfy all the constraints 
end if 

if no non-dominated solutions can be found 
Generate random weights Wk {k = 1, ...,K) 

Find the best solution g\ (among all fireflies) to minimize in (4) 
Random walk around using (5) 
end if 

Update and pass the non-dominated solutions to next iterations 
end 

Sort and find the current best approximation to the Pareto front 
Update t ^ t-\-l 
end while 

Postprocess results and visualisation; 

Figure 1: Pseudo Code: Multiobjective firefly algorithm (MOFA). 

The procedure starts with an appropriate definition of objective functions with associated non- 
linear constraints. We first initialize a population of n fireflies so that they should distribute among 
the search space as uniformly as possible. This can be achieved by using sampling techniques via 
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uniform distributions. Once the tolerance or a fixed number of iterations is defined, the iterations 
start with the evaluation of brightness or objective values of all the fireflies and compare each pair 
of fireflies. Then, a random weight vector is generated (with the sum equal to 1), so that a com- 
bined best solution g\ can be obtained. The non-dominated solutions are then passed onto the next 
iteration. At the end of a fixed number of iterations, in general n non-dominated solution points 
can be obtained to approximate the true Pareto front. 

In order to do random walks more efficiently, we can find the current best which minimizes a 
combined objective via the weighted sum 

K K 

"^i^) = X! ^kfk, ^ Wfe = 1. (4) 
fe=l fc=l 

Here Wk = Pk/K where pk are the random numbers drawn from a uniform distributed Unif[0,l]. In 
order to ensure that X^j.«^/c = 1, a rescaling operation is performed after generating K uniformly 
distribiitcd numbers. It is worth pointing out that the weights Wk should be chosen randomly at 
each iteration, so that the non-dominated solution can sample diversely along the Pareto front. 
If a firefly is not dominated by others in the sense of Pareto front, the firefly moves 

xl+^ =gl + aiel (5) 

where g* is the best solution found so far for a given set of random weights. 

Furthermore, the randomness can be reduced as the iterations proceed, and this can be achieved 
in a similar manner as that for simulated annealing and other random reduction techniques [11]. We 
will use 

at = Q!o0.9*, (6) 

where ao is the initial randomness factor. 



2.3 Pareto Optimal Front 

For a minimization problem, a solution vector u = (ui, ..,Un)'^ is said to dominate another vector 
V = {vi, ...,Vn)'^ if and only if Ui < Vi for \/i S {1, ...,n} and 3i G {1, ...,n} : Ui < Vi. In other words, 
no component of u is larger than the corresponding component of v, and at least one component is 
smaller. Similarly, we can deflne another dominance relationship ^ by 

u ^ V u -< vW u = V. (7) 

It is worth pointing out that for maximization problems, the dominance can be defined by replacing 
-< with Therefore, a point is called a non-dominated solution if no solution can be found that 
dominates it [5]. 

The Pareto front PF of a multiobjective can be defined as the set of non-dominated solutions 
so that 

PF = {seS^s' eS:s' ^s}, (8) 

where <S is the solution set. 

To obtain a good approximation of the Pareto front, a diverse range of solutions should be 
generated using efficient techniques [15, 29, 30, 31]. For example. Levy flights ensure the good 
diversity of the solutions, as we can see from later simulations. 



3 Numerical Results 

We have implemented the proposed MOFA in Matlab, and we have first validated it against a sot of 
multiobjective test functions. Then, we have used it to solve some industrial design of structures. In 
order to obtain the right algorithm-dependent parameters, we have carried out detailed parametric 
studies. 
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3.1 Parametric Studies 



By varying the parameters ao, a-nd 7, we have carried out parametric studies by setting ao = 
to 1 with a step of 0.05, /3o = to 1 with a step of 0.05, and 7 = 0.1 to 10 with a step of 0.1 and 
then 1. From simulations, we found that we can use ao = 0.1 to 0.5, (3o = 0.7 to 1.0 and 7 = 1 for 
most problems. 

The stopping criterion can be defined in many ways. We can either use a given tolerance or a 
fixed number of iterations. For a given tolerance, we should have some prior knowledge of the true 
optimum of the objective function so that we can calculate the differences between the current best 
solutions and the true optimal solution so that we can assess if the tolerance is met. In reality, we 
usually do not know the true optimum in advance, except for a few well-tested cases. In addition, 
the number of functions can vary significantly from function to function even for the same tolerance. 
From the implementation point of view, a fixed number of iterations is not only easy to implement, 
but also suitable to compare the closeness of Pareto front of different functions. So we have set 
the fixed number iterations as 2500, which is sufficient for most problems. If necessary, we can also 
increase it to a larger number. 

One can generate points of the Pareto front in two ways: increase the population size n or run 
the program a few more times. Through simulations, we found that increasing n typically leads to a 
longer computing time than re-running the program a few times. This may be due to the fact that 
manipulations of large matrices or longer vectors usually take longer. In order to generate M points 
using a smaller population size n, it requires to run the program M/n times, each run with different, 
random initial configurations but with the same number of iterations t = 2500. For example, to 
generate M = 200 points, we can use n = 50, which is easily done within a few minutes. Therefore, 
in all our simulations, we will use the fixed parameters: n = 50, ao = 0.25, /3o = 1 and 7 = 1. 

3.2 Multiobjective Test Functions 

There are many different test functions for multiobjective optimization [32, 33, 34], but a subset of 
a few widely used functions provides a wide range of diverse properties in terms of Pareto front and 
Pareto optimal set. To validate the proposed MOFA, we have selected a subset of these functions 
with convex, non-convex and discontinuous Pareto fronts. We also include functions with more 
complex Pareto sets. To be more specific in this paper, we have tested the following 5 functions: 

• Schaffer's Min-Min (SCH) test function with convex Pareto front [21, 35] 

h{x) = x^, f2{x) = {x-2f, -103<a;<103. (9) 

• ZDTl function with a convex front [33, 34] 



fi{x)=xi, f2{x) = g{l- Vh/g), 

9=^ + ^fZY^^ X, e [0,1], z = l,...,30, (10) 

where d is the number of dimensions. The Pareto optimality is reached when 5 = 1, and thus 

the true Pareto front is /2 = 1 — V/i • 

ZDT2 function with a non-convex front 

ZDT3 function with a discontinuous front 

fi{x) = xi, f2{x) = g 1- ~ sin(107r/i) 

where g in functions ZDT2 and ZDT3 is the same as in function ZDTl. In the ZDT3 function, 
/i varies from to 0.852 and /2 from -0.773 to 1. 
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• LZ function [20, 36] 

fi= xi + — — ^ \xj - sin(67ra;i + ^) , 

f2 = l-V^+ +^ k" ~ sin(67ra;i + , (11) 

|J2| ^.g^^ a 

where Ji = is odd } and J2 = {j\j is even } where 2 < j < d. This function has a Pareto 
front /2 = 1 — -s/Ti with a Pareto set 

= sin(67ra;i + ^), j = 2,3, a;ie[0,l]. (12) 

After generating 200 Pareto points by MOFA, these points are compared with the true front /2 = 
1 - ^/h of ZDTl (see Fig. 2). 

Let us define the distance or error between the estimated Pareto front PF"^ to its corresponding 
true front PF* as 

JV 

Ef = WPF"" - PF*\\^ = ^(PF/ - p^;*)^ (13) 

where N is the number of points. For all the test functions, the true Pareto fronts have analyti- 
cal forms [20, 33, 35], for example, /2 = 1 ^ vTi for ZDTl, which makes the above calculations 
straightforward. 

The convergence property can be viewed by following the iterations. As this measure is an 
absolute measure, which depends on the number of points. Sometimes, it is easier to use relative 
measure using generalized distance 

N 

Y,{PF,-PFjY. (14) 

Fig. 3 shows the exponential-like decrease of Dg as the iterations proceed. We can see clearly 
that our MOFA algorithm indeed converges almost exponentially. The results for all the functions 
are summarized in Table 1, and the estimated Pareto fronts and true fronts of other functions are 
shown in Fig. 4 and Fig. 5. In all these figures, the vertical axis is /2 and the horizontal axis is /i. 



7 




200 400 600 800 1000 

iterations 

Figure 3: Convergence of the proposed MOFA. The least-square distance (vertical axis) from the 
estimated front to the true front of ZDTl for the first 1000 iterations. 



Table 1: Summary of results. 



Functions 


Errors (1000 iterations) 


Errors (2500 iterations) 


SCH 


5.5E-09 


4.0E-22 


ZDTl 


2.3E-6 


5.4E-19 


ZDT2 


8.9E-6 


1.7E-14 


ZDT3 


3.7E-5 


2.5E-11 


LZ 


2.0E-6 


7.7E-12 



3.3 Comparison Study 

In order to compare the performance of the proposed MOFA with other established multiobjec- 
tive algorithms, we have carefully selected a few algorithms with available results from the lit- 
erature. When the results have not been available, we have implemented the algorithms using 
well-documented studies and then generated new results using these algorithms. In particular, we 
have used other methods for comparison, including vector evaluated genetic algorithm (VEGA) 
[21], NSGA-II [37], muhiobjective differential evolution (MODE) [23, 38], differential evolution for 
multiobjective optimization (DEMO) [22], muhiobjective bees algorithms (Bees) [39], and strength 
Pareto evolutionary algorithm (SPEA) [37, 40]. The performance measures in terms of generalized 
distance Dg are summarized in Table 2 for all the above major methods. 

It is clearly seen from Table 2 that the proposed MOFA obtained better results for almost all 
five cases, though for ZDT2 function our result is the same order (still slightly better) as that by 
DEMO. 

4 Design Optimization 

Design optimization, especially design of structures, has many applications in engineering and in- 
dustry. As a result, there are many different benchmarks with detailed studies in the literature 
[39, 41, 42]. Some benchmarks have been solved by various methods, while others do not have all 
available data for comparison. Thus, we have chosen the welded beam design, and disc brake design 
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Figure 4: a) Pareto front of test function SCH, and b) Pareto front of test function ZDT2. 




Figure 5: a) Pareto front of test function ZDT2, and b) Pareto front of test function LZ. 



Table 2: Comparison of Dg for n = 50 and t = 500 iterations. 



Methods 


ZDTl 


ZDT2 


ZDT3 


SCH 


LZ 


VEGA 


3.79E-02 


2.37E-03 


3.29E-01 


6.98E-02 


1.47E-03 


NSGA-II 


3.33E-02 


7.24E-02 


1.14E-01 


5.73E-03 


2.77E-02 


MODE 


5.80E-03 


5.50E-03 


2.15E-02 


9.32E-04 


3.19E-03 


DEMO 


1.08E-03 


7.55E-04 


1.18E-03 


1.79E-04 


1.40E-03 


Bees 


2.40E-02 


1.69E-02 


1.91E-01 


1.25E-02 


1.88E-02 


SPEA 


1.78E-03 


1.34E-03 


4.75E-02 


5.17E-03 


1.92E-03 


MOFA 


1.90E-04 


1.52E-04 


1.97E-04 


4.55E-06 


8.70E-04 
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among the well-known benchmarks [42, 43, 44]. In the rest of this paper, we will solve these two 
design benchmarks using MOFA. 



4.1 Welded Beam Design 



Multiobjective design of a welded beam is a classical benchmark which has been solved by many 
researchers [6, 12, 44]. The problem has four design variables: the width w and length L of the 
welded area, the depth d and thickness h of the main beam. The objective is to minimize both the 
overall fabrication cost and the end deflection S. 

The detailed formulation can be found in [6, 12, 44]. Here we only rewrite the main problem as 



minimise fi{x) = 1.10471^^^ + 0.04811d/i(14.0 + L), minimize /a = 6, 



(15) 



subject to 



where 



9i{x 
92{x 
93{x 

9g{x 
97{x 



= w-h<0, 

= 6{x) - 0.25 < 0, 

= t{x) - 13, 600 < 0, 

= <7{x) - 30,000 < 0, 

= 0.10471W72 + 0.04811M(14+ L) 

= 0.125- w<0, 

= 6000 - P{x) < 0, 



(16) 



5.0 < 0, 



a{x) = 



_ 504,000 



D = l^L^ + {w + d)^, 



g = 6000(14 +f), 



65,856 
30,000/td3 ' 



QD 
J ' 



(17) 



6000 
V2wL' 



t{x) 



P = 0.61423 X 10*^ ^ 



V-*- 28 /• 



The simple limits or bounds are 0.1 < i, d < 10 and 0.125 < w,h < 2.0. 

By using the MOFA, we have solved this design problem. The approximate Pareto front generated 
by the 50 non-dominated solutions after 1000 iterations are shown in Fig. 6. This is consistent with 
the results obtained by others [39, 44]. In addition, our results are more smooth with fewer iterations, 
which shows the efficiency of the proposed MOFA. A comparison of the results with those obtained 
by other methods is shown in Fig. 7 where we can see MOFA converged faster. 



4.2 Design of a Disc Brake 

Design of a multiple disc brake is another benchmark for multiobjective optimization [12, 18, 44]. 
The objectives are to minimize the overall mass and the braking time by choosing optimal design 
variables: the inner radius r, outer radius R of the discs, the engaging force F and the number of the 
friction surfaces s. This is under the design constraints such as the torque, pressure, temperature, 
and length of the brake. This bi-objective design problem can be written as: 

Minimize Mx) = 4.9 x IQ-^R' - r'){s - 1), f2{x) = %^(j^3 _ ^3) ' (18) 
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Figure 6: Pareto front for the bi-objective beam design where the horizontal axis corresponds to 
cost and the vertical axis corresponds to deflection. 
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subject to 

gi{x) = 20- (R-r) <0, 
.g2(a;) = 2.5(s + l)-30<0, 

.93(3;) -^Tllfc^y- 0-4 <0, 

19 

g,{x) = - 1 < 0, 

35(0^) = 900 -MWM<o. 

The simple limits are 

55 < r < 80, 75 < i? < 110, 1000 <F < 3000, 2 < s < 20. (20) 

The detail formulation of this problem and background descriptions can be found in [41, 12, 18, 44]. 
The Pareto front of 50 solution points after 1000 iterations obtained by MOFA is shown in Fig. 8, 
where we can see that the results are smooth and are the same or better than the results obtained 
in [44]. This can also be seen from the comparison of converge rates shown in Fig. 9. 

soF. 

25 . 

20 



CM 



15 
10 



■ 



Figure 8: Pareto front for the disc brake design where /2 is the vertical axis and /i is the horizontal 
axis. 

In order to see how the proposed MOFA performs for the real-world design problems, we also 
solved the same problems using other available multiobjective algorithms. The comparison of the 
convergence rates is plotted in the logarithmic scales in Fig. 7 and Fig. 9. We can see from these 
two figures that the convergence rates of MOFA are of the highest in an exponentially decreasing 
way in both cases. This again suggests that MOFA provides better solutions in a more efficient way. 

The simulations for these benchmarks and test functions suggest that MOFA is a very efficient 
algorithm for multiobjective optimization. It can deal with highly nonlinear problems with complex 
constraints and diverse Pareto optimal sets. 



5 Conclusions 

Multiobjective optimization problems are typically very difficult to solve. We have successfully 
formulated a new algorithm for multiobjective optimization, namely, multiobjective firefly algorithm 
(MOFA), based on the recently developed firefly algorithm. The proposed MOFA has been tested 
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Figure 9: Convergence comparison for the disc brake design. 

against a subset of well-chosen test functions, and then been applied to solve design optimization 
benchmarks in industrial engineering. 

By comparing with other algorithms, the present results suggest that MOFA is an efficient 
multiobjective optimizer. Further studies can focus on parametric studies which can be very useful 
to identify the optimal ranges of parameters for various optimization problems. 

In addition, convergence analysis of MOFA will also show insight in the working mechanism of 
the algorithm, and this may also help to improve the proposed algorithm or even design new algo- 
rithms. For example, hybridization with other algorithms may also prove to be fruitful. Furthermore, 
formulation of a discrete MOFA will also be an important topic for further research. 
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